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ABSTRACT 

We present a near-IR spectrum of the nearby radio galaxy PKS1549-79 (z = 0.153). 
These data were taken with the aim of testing the idea that this object contains a 
quasar nucleus that is moderately extinguished, despite evidence that its radio jet 
points close to our line-of-sight. We detect broad Paa emission (FWHM 1745 ± 40 
km s _1 ), relatively bright continuum emission, and a continuum slope consistent with 
a reddened quasar spectrum (3.1 < Ay < 7.3), all emitted by an unresolved point 
source. Therefore we conclude that we have, indeed, detected a hidden quasar nucleus 
in PKS1549-79. Combined with previous results, these observations are consistent 
with the idea that PKS1549-79 is a young radio source in which the cocoon of debris 
left over from the triggering events has not yet been swept aside by circumnuclear 
outflows. 
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1 INTRODUCTION 

In recent years, much AGN research has been con ducted in 
the context of "unified schemes" fe.g. lBarthell989T) : i.e. that 
differences between certain types of active galaxy can be ex- 
plained by orientation effects. A great deal of AGN research 
is based around refining and constraining such models. How- 
ever, the simplest unified models generally describe a static 
and unchanging situation, a steady-state regime, and this is 
unlikely to be realistic. As the activity evolves it is likely 
to affect the distribution of ISM surrounding the nucleus, 
potentially influencing the evolution of the host galaxy. To 
better understand these processes it is necessary to develop 
a more dynamic, evolutionary model for AGN, and to this 
end it is important to study sources at different evolutionary 
stages. 

The existence of compact steep-spectrum (CSS) and 
gigahertz-peaked (GPS) radio sources, which are intrinsi- 
cally small, illustrates the fact that there is a large range in 
the physical dimensions of extragalactic radio sources. The 
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compactness of these sources must either be due to frus- 
tration, where a dense ISM is inhibiting the expansion of 
the radio lobes, or to youth, in which case the radio lobes 
simply have not yet had time to expand. Of these two poss i- 
bilities youth seems the most likely (e.g. lFanti et al.l fl995). 
Broad forbidden lines are also much more common in these 
compact sources, sugge sting that the cores of th ese sources 
are more un s ettled dGelderman fc Whittlell994t) . Moreover, 
iBaker et al.l J2002) find an anti-correlation between C IV 
absorption - which they show to be associated with dusty 
regions - and radio source dimensions; the absorption be- 
coming less pronounced with increasing size. It is therefore 
becoming apparent, and perhaps expectedly so, that young 
radio-loud AGN suffer greater dust extinction than their 
more evolved counterparts, with the nucleus still surrounded 
by a cocoon of dust and gas. 



Based on its emission line widths and kinematics, the 
southern radio galaxy PKS1549-79 (z = 0.153) appears to 
be just such a young source. The observed properties of this 
unusu al object are discussed in detail in iTadhunter et alJ 
(2001) (hereafter T2001). The flat spectrum, compactness 
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of the radio emission (~150 milliarcsec, or ~540 pc ) and 
one-sided jet morphology of this source indicate that its ra- 
dio axis is aligned close to our line-of-sight. In the standard 
unification model, the source would be expected to resem- 
ble a quasar at optical wavelengths, with non-stellar contin- 
uum emission and broad permitted lines. However, optical 
spectra of PKS1549-79 show a predominantly stellar contin- 
uum with no good evidence for broad permitted lines, but 
with strong and unusually broad forbidden lines (~1300 km 
s _1 ). Amongst powerful radio galaxies such broad forbidden 
lines are only ever detected in sources that are intrinsically 
compact. The blueshift of the [O m]AA5007, 4959 emission 
lines relative to low-ionisation forbidden lines and the Hi 21 
cm absorption line suggests that the [O ill] emitting region 
is ou tflowing at 600 ± 60 km s" 1 (T2001, iMorganti et alJ 
1200 lib 

Overall, observations suggest that there is significant 
dust obscuration in this object even though our line-of-sight 
is close to the radio axis. This is entirely consistent with 
PKS1549-79 being a young radio source in which the ob- 
scuring material around the nucleus has not yet been swept 
aside (T2001). If this model is correct, near-infrared observa- 
tions - being subject to only ~10% of the visual extinction - 
should show a bright, non-stellar continuum and broad per- 
mitted lines as the quasar shines through the ISM. This is 
the prediction we test in this paper. 



2 DATA COLLECTION AND REDUCTION 

The K-band infrared spectra were taken in shared risks ser- 
vice mode on 26 July 2002 using the IRIS2 instrument on 
the AAT. The Sapphire 240 grism was used with a 1 arc- 
sec slit, oriented North-South on the sky. The seeing was 
reported to be 1.1 arcsec throughout the observations; anal- 
ysis of stars in the 10 second exposure acquisition image is in 
good agreement with this (see Section 3.1). 22 exposures of 
300 seconds were obtained for PKS1549-79; the galaxy was 
'nodded' between two positions on the slit in an ABBA pat- 
tern with 11 exposures at each position. Four exposures of 
80 seconds were taken of the AO star HIP77712 at a similar 
airmass for calibration purposes. 

Each set of 11 galaxy exposures was co-added in the 
iraf package using median filtering to remove cosmic rays; 
the median 'B' image was then subtracted from the me- 
dian 'A' image to remove the night sky lines. The standard 
star frames were combined in a similar way. The resulting 
galaxy and star frames were then flat-fielded using a dark- 
subtracted flat field taken with the same instrumental setup. 

A xenon arc frame was used to make a two-dimensional 
wavelength calibration which was then applied to the data 
frames. The original spectra cover the range 20109 - 24640 
A but the useful data presented here are over the range 
20500 - 23000 A. From the night sky lines, the accuracy of 
the wavelength calibration was found to be ± 0.65 A (± 
9 km s _1 ) and the systematic errors are estimated to be < 
1 A (< 14 km s _1 ). Fits to the night sky lines show the 
spectral resolution to be 10.1 ± 0.4 A; the pixel scale of 4.43 

1 We assume cosmological parameters of Hq = 50 km s -1 Mpc - 1 
and qo = throughout. 
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Figure 1. K-band spectrum of PKS1549-79 obtained with the 
IRIS2 instrument on the AAT. Doubling of the lines is caused by 
birefringence associated with the grism used for the observations 
(See Section 2). 

A pix -1 is adequate to sample this. The spatial pixel scale of 
the acquisition image is 0.446 arcsec pix -1 . One-dimensional 
spectra were extracted from the two-dimensional frames, 
each from 11 pixel rows to aid accurate flux calibration. 

The exposures of HIP77712 were used to flux calibrate 
the data, with the assumption that the intrinsic spectral en- 
ergy distribution (SED) of the star is that of a perfect black 
body at T = 9480 K. The m agnitude-to-flux convers i on wa s 
performed with reference to lBessell. Castelli fc Plea ([1998). 

During data reduction it became apparent that there 
was a problem with the IRIS2 instrument in spectroscopy 
mode at the time of the observations. All the spectra we had 
obtained, including arc calibration frames, showed a split- 
ting of spectral features in the wavelength direction. After 
consulting with the AAT team we were informed that the 
manufacturers had mounted the Sapphire grisms at right- 
angles to their intended orientation and that a birefringence 
effect was leading to the observed doubling. The spectra 
we obtained are composites of the actual spectrum and a 
copy of this spectrum, blueshifted by around 50 pixels (~230 
A in the case of these data) relative to the original and at a 
slightly lower intensity. However, the splitting caused by the 
birefringence is large enough that useful kinematic informa- 
tion can be extracted for individual lines in the spectra. The 
final, reduced K-band spectrum of PKS1549-79 is shown in 
Fig. The Starlink FIGARO and DIPSO packages were used 
to analyse the data. 



3 RESULTS 

3.1 Paschen alpha 

A single Gaussian provides an adequate fit to each of the 
two components of the split line. The best fit to the Paa 
line(s) is shown in Fig. |5J with the redder peak being at the 
correct wavelength. A constraint was applied that the two 
Gaussians should have the same width. The results of the 
fitting are summarised in Tabic 

The fit gives a FWHM of 126 ± 3 A (1745 ± 40 
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Parameter Pact H2 or [Si vi] 

FWHM (km s" 1 ) 1745 ± 40 1100 ± 180 

z 0.15266 ± 0.00007 0.1525 ± 0.0002 or 0.1494 ± 0.0002 

Total line flux (erg cm" 2 s" 1 ) (3.16 ± 0.08) X 10" 14 (3.5 ± 1.0) X 10" 15 

Continuum flux (erg cm" 2 s" 1 A" 1 ) (4.0 ± 0.1) X 10" 16 (4.0 ± 0.2) xl0~ 16 

Equivalent width (A) 79 ± 3 9 ± 3 

Rest-frame equivalent width (A) 69 ± 3 8 ± 3 



Table X. Results of the DIPSO fits to the spectral lines. The instrumental width has a negligible effect on the measured FWHMs. 
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Figure 2. DIPSO fit to the Pact line(s). Pace wavelengths corre- 
sponding to the [O 11] and [O in] redshifts (T2001) are marked 
for comparison. 



km s 1 in the rest frame). This value is consistent with 
the widths of permitted li nes in broad line radio galax- 
ies (BLRG) and qu asars jHill. Goodrich fc DePovl Il99fl 
iMarziani et al-lfeoOll) . although it is at the lower end of the 
range. The relatively small broad-line widths are also con- 
sistent with_t^ie_£adjoiet_being oriented close to our line-of- 
sight JWills fc Brownelll98d) . 

The redshift of the Pact line is found to be z = 0.15266 
± 0.00007. This is consistent with the [O 11] redshift of z = 
0.1526 ± 0.0002 but inconsistent with the [O in] (outflow) 
redshift of z = 0.1501 ± 0.0002 (T2001). This indicates that 
the Paa emission originates in a broad line region (BLR) 
at the systemic redshift. The Paa wavelengths correspond- 
ing to the [O 11] and [O in] redshifts are marked on |5] for 
comparison. 

The K-band continuum flux level is ~4x that measured 
in the optical spectrum of T2001, despite the fact that the 
optical spectrum was obtained using a larger aperture (4.3 
by 5 arcsec). This is inconsistent with the SED of an unred- 
dened stellar population, particularly a young population, 
but is consistent with that of a moderately obscured quasar 
source shining through the intervening ISM. A power-law 
was fit to the continuum to aid comparisons with quasar 
continua (see Section 4). The spectral index of this power- 
law is a = 2.3 ± 0.1 (F„ tx v~ a , F A oc A a ~ 2 ). This index 
is significantly larger than that measured f or the 16 unred- 
dened sources in the lSimpson fc Rawlingsl i2000T) sample of 
radio- loud quasars (-0.67 < a < 1.62), but similar to the two 



reddened quasars in that sample. If the intrinsic spectrum 
of PKS1549-79 is quasar-like then this indicates significant 
reddening. 

The combined flux of both Paa components is (3.16 
± 0.08) x 10~ 14 erg cm -2 s _1 . This flux gives a rest-frame 
equivalent width of 69 ± 3 A for an assumed continuum flux 
of (4.0 ± 0.1) x 10~ 16 erg cm -2 s _1 A -1 in the observed 
frame. For comparison, the rest-frame Paa equivalent widths 
measured in the low-redshift quasars 3C 2 73 and PDS 456 
are 167 ± 16 A and 120 ± 8 A respectively (iHill et alJll996l. 
ISimpson et alll999l) . Both of these quasar values are signifi- 
cantly larger than that measured in PKS1549-79. This could 
indicate the presence of a beamed, non-thermal continuum 
component in PKS1549-79 which would further support the 
theory that the radio jet is oriented towards our line-of-sight. 

If there is a quasar nucleus shining through at near- 
IR wavelengths then the emission from the galaxy should 
be dominated by a central, unresolved source. To test this, 
the spatial intensity profile of PKS 1549-79 was compared 
with those of four bright stars in the acquisition image. The 
mean FWHM of the stellar 2-D Gaussian measures 1.13 ± 
0.06 arcsec in the x-direction by 1.23 ± 0.06 arcsec in the 
y-direction. In comparison, the profile of PKS1549-79 mea- 
sures 1.13 ± 0.17 by 1.43 ± 0.15 arcsec. Therefore PKS1549- 
79 appears quasar-like at near-IR wavelengths. 



3.2 Other features 

Apart from the Paa line, the only other feature in the 
spectrum is a second line visible at ~22550 A, subject to 
the same splitting as the rest of the spectrum. As the line 
appears faint, the spectrum was re-binned to one-quarter 
resolution over this region to improve the signal to noise. 
The line is broad enough for this to be acceptable. The 
blueshifted split component is too noisy for an adequate fit 
to be made and so only the redmost component was fit in 
DIPSO. The fit is shown in Fig. Eland the results of the fitting 
are given in Table 

The FWHM of this line is 83 ± 14 A (1100 ± 180 km 
s _1 in the rest frame). Correcting for the fact that only 
one component was fit, the total flux of the line is (3.5 ± 
1.0) x 10 -15 erg cm -2 s _1 . This corresponds to an isotropic 
luminosity of (4.1 ± 1.2) x 10 41 erg s" 1 . 

The most likely candidate for this emission is the H2 v 
= 1-0 S(3) line, our identification is based on redshift argu- 
ments. The H2 v = 1-0 line would lie at the same redshift 
as the Paa line, whereas other lines imply large inflow or 
outflow velocities. However, [Si vi] is also a possibility since 
the implied outflow velocity and width are broadly consis- 
tent with the [O ill] outflow parameters (T2001). The pos- 
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Figure 3. DIPSO fit to the -22550 A line. 



sibility that this feature is a blend of both lines cannot be 
discounted. However, the low S/N of the data - coupled with 
the uncertainties introduced by the splitting of the spectrum 
- means that no definite conclusions can be drawn. 

If the feature is indeed H 2 v = 1-0 S(3), then PKS1549- 
79 is one of the most luminous molecular hydrogen emit- 
ters in the local universe, surpassing the H2 l uminosity of 
Cygnus A by almost an order of magnitu de (|Ward et alJ 
ll99ll.lThornton. Stockton fc Ridgwavlll999D . 

Unfortunately no other molecular hydrogen lines are 
visible in the spectrum with which to compare this result. 
This is to be expected however, as the H2 v = 1-0 S(3)/[Si 
vi] feature is only weakly detected and the only molecular 
line expected to be comparable in luminosity (H2 v = 1-0 
S(l)) is redshifted out of the useful range of our data. 



4 DISCUSSION 

Our observations provide clear evidence that emission from 
the central regions of PKS 1549-79 is significantly affected by 
dust obscuration. In this section we attempt to estimate the 
degree of reddening 2 and correct for it in order to obtain a 
more realistic estimate of the intrinsic Paa luminosity. For 
comparison, T2001 derive a visual extinction in the range 
0.23 < Ay < 18 mag, based on Hi 21 cm absorption char- 
acteristics. 

Firstly, we try to determine the reddening directly from 
Balmer line ratios, these being relatively insensitive to vari- 
ations in the phy sical conditions of the emitting medium. 
With reference to lTadhunter et alJ lll993l) . the Paa/H/3 line 
ratio is found to be 25 ± 1. Assuming Case B recombina- 
tion, a temperature of 10,000K and an electron density of 
10 4 e lectrons cm" 3 , this line ratio g ives an Es_v = 1.29 ± 
0.01 iSeatonlll979t IOsterbrocklll989r) . corresponding to a vi- 
sual extinction of Ay = 4.00 ± 0.13 mag (Av = R x Eb-v, 
assuming R = 3.1 ± 0.1). However, there is a great deal of 
uncertainty associated with this value because the H/3 emis- 
sion is likely to be dominated by a narrow line region (NLR) 



component. The BLR emission from within the obscuring re- 
gion may then contribute only a fraction, if any, of the total 
H/3 flux, and therefore the derived extinction is likely to be 
an underestimate. Recently obtained high-resolution spectra 
of PKS1549-79 include Ha emission: however, the Ha line 
is contaminated with bright N n recombination lines and an 
atmospheric absorption line. These features will significantly 
reduce the accuracy of line fittings to the Ha line and for 
this reason the line has not been used in this analysis. 

Another approach is to fit a power-law to the near-IR 
continuum and then calculate the de-reddening necessary to 
make the slope consistent with a quasar spectrum. As stated 
in Section 3.1, the spectral index of the power-law fit to our 
data is a = 2.3 ±0.1. Typical unreddened quasar spectral 
indices lie in the range -0.67 < a < 1.62 and average a 
= 0.90 dSimpson fc Rawlingsl I2OO0T) . Clearly the spectrum 
for PKS1549-79 is considerably redder than even the upper 
limit. De-reddening of the spectrum to match the quasar 
range of values requires a visual extinction in the range 1.8 
< Av < 8.0 with Av = 3.8 for the average case. The full 
de-reddening results are given in Table If the spectrum 
is contaminated by a beamed, non-thermal component then 
the intrinsic continuum shape may not resemble a quasar; 
in the extreme case it may more closely rese mble a BL Lac. 
The 6 BL Lacs in the iMassaro et alJ 1119951) sample show a 
tighter range of spectral indices (0.86 < a < 1.19) with an 
average that is redder than that of the quasar sample: this 
may imply a smaller degree of reddening. However, as the 
full range of BL Lac spectral indices lies within the broader 
range of quasar spectral indices, this result has no significant 
effect on our analysis. 

To investigate the contribution such reddened power- 
law continua would make at optical wavelengths, curves were 
generated for the derived power-laws and then re-reddened 
by the corresponding amount. Fig.0]shows the results of this 
for the maximum, average and minimum extinctions/power- 
laws derived above. Certainly, the minimum reddening case 
is ruled out on the basis that the quasar alone would con- 
tribute more flux at >6000 A than is detected (T2001). The 
average extinction, corresponding closely with that derived 
from the Paa/H/3 ratio, would also be contributing around 
half the observed flux in this region and also with a marked 
slope that one would expect to detect superimposed on the 
starlight: modelling sh ows no evidence for a p ower-law com- 
ponent in this region jTadhunter et al . 2002). This implies 
that the actual reddening is closer to the upper limit than 
to the average, regardless of whether the intrinsic contin- 
uum slope more closely resembles a quasar or BL Lac. It 
seems safe to conclude that the total visual extinction lies 
in the range 3.8 < Ay < 8.0 mag. This corresponds to an 
intrinsic extinction in the host galaxy of 3.1 < Av < 7.3, 
having corrected for the Galactic e xtinction of Ay = 0.68 
JSchleeel. Finkbeiner fc DavisllT998T) . 

The de-reddened Paa fluxes for the three cases are given 
in Table |5J along with the corresponding isotropic lumi- 
nosities. Comparison with similarly derived luminosities for 



2 We assume the standard interstellar extinction law of W hitfordl 
(1958). 



a sample of qu asars and BLRGs llRudv fc Tokunagalll98l 
iHill et al]|l996l and lSimpson" ct al. 1999) shows that the lu- 
minosity corresponding to the average reddening ((59.0 ± 
1.5) x 10 41 erg s _1 ) is more than an order of magnitude less 
than that measured in the low-z quasars, and at the upper 
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Reddening Degree > Minimum Average Maximum 

Spectral Index a 1.62 0.90 -0.67 

E B _ V 0.60 1.23 2.62 

A v (total) 1.8 3.8 8.0 

A v (intrinsic)" 1.2 3.1 7.3 

Pact Flux (10~ 14 erg s" 1 A" 1 ) 3.95 ± 0.10 5.02 ± 0.13 8.50 ± 0.22 

Pact Luminosity 6 (10 41 erg s" 1 ) 46.4 ± 1.2 59.0 ± 1.5 99.8 ± 2.5 



Table 2. Resulting Paa emission parameters corresponding to the different degrees of reddening as discussed in the text. 
"Assuming Galactic extinction Ay = 0.68 ISchlegel et alJl99Sft . 
^Calculated assuming isotropic emission. 




Wavelength (A) 

Figure 4. Reddened power-laws generated for the parameters 
given in Table [5] The average case is shown by the solid line, 
the minimum case by the upper broken line and the maximum 
case by the lower broken line. The thick solid lines represent the 
optical continuum data from T2001 and the near-IR data from 
this paper. The data are plotted at the observed wavelengths and 
the model shifted into this frame to match. 

end of the range measured in the BLRGs. These results are 
consistent with PKS1549-79 being a low-luminosity quasar. 



5 CONCLUSIONS AND FURTHER WORK 

The detection of a broad Paa line and a bright, reddened 
infrared continuum indicates that we have, indeed, uncov- 
ered an obscured BLR and quasar nucleus in PKS1549-79. 
The calculated intrinsic luminosity and equivalent width 
of the Pact emission, and the star-like profile of the near- 
IR nucleus, further support this co nclusion. Although such 
findings are n ot unprecedented llDiorgovskv et alJ Il99ll 
iHill et aljll99r3) . PKS1549-79 is a special case since the ob- 
served radio properties lead us to believe that the AGN 
is viewed from a direction close to that of the radio axis 
and should therefore exhibit the optical/infrared properties 
of radio loud quasars; i.e. no significant obscuration of the 
BLR and nucleus. The fact that significant obscuration is 
seen suggests that PKS1549-79 is a fundamentally different 
type of object that does not fit neatly into standard unified 
schemes. 

It is our suggestion that PKS 1549-79 is an AGN in the 



early stages of evolution, a proto-quasar. The high obscu- 
ration is a transitory phase that will pass as the gas and 
dust is dissipated or blown out of the ionisation cones by 
the circumnuclear outflows detected at optical wavelengths 
(T2001). These infrared data support this model and tie in 
well with expectations. 

Since a key assumption of the model is that the radio 
jet is oriented towards us, it is important to validate this 
assumption. The fact that the Paa equivalent width is no- 
ticeably smaller than that for a typical quasar is consistent 
with the presence of a beamed continuum component and 
may therefore support this assumption. However, this is not 
sufficient evidence in itself. Multi-epoch VLBI imaging of 
the radio emission could provide incontrovertible evidence 
that the radio axis is oriented towards us by revealing super- 
luminal motion in the jet. 
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